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Abstract 

Dynamical system methods are used in the study of the stabihty of spatially flat homogeneous 
cosmologies within a large class of generalized modified gravity models in the presence of a relativistic 
matter-radiation fluid. The present approach can be considered as the generalization of previous works 
in which only F{R) cases were considered. Models described by an arbitrary function of all possible 
. geometric invariants are investigated and general equations giving all critical points are derived. 

o 

. 1 Introduction 
D . 

I . Recent cosmological data support the fact that there is a strong evidence for a late accelerated expansion 

of the observable universe, apparently due to the presence of an effective positive and small cosmological 
constant of unknown origin. This is known as dark energy issue (see for example pLJ). 

Modified gravity models are possible realizations of dark energy (for a recent review and alternative 
approaches see [HH]), which may offer a quite natural geometrical approach again in the spirit of the 
original Einstein theory of gravitation. In fact, the main idea underlying these approaches to dark 
, energy puzzle is quite simple and consists in adding to the gravitational Einstein-Hilbcrt action other 

^ ' gravitational terms which may dominate the cosmological evolution during the very early or the very late 

, universe epochs, but in such a way that General Relativity remains valid at intermediate epochs and also 

at non cosmological scales. They are generalization of the ACDM model, namely Einstein gravity plus a 
positive cosmological constant, the simplest model for dark energy, which however has to be confronted 
^ ' with well known difficulties, among them, the cosmological constant problem, an unsolved issue so far. 

i In the present paper, first of all we shall generalise the analysis presented in A.^, where the stability 

i of the de Sitter solution (the vacuum invariant submanifold) has been investigated in a class of modified 

gravitational cosmological models defined in a Friedmann-Robertson- Walker (FRW) spatially flat space- 
time. The method is well known and consists in rewriting the equations of motion as a system of first-order 
I autonomous differential equations and makes use of the theory of dynamical systems (see El H |H1 E] 

and references therein). We remind that the stability or instability issue is really relevant in cosmology. 
I For example, in the ACDM model it ensures that no future singularities will be present in the solution. 

Within cosmological models, the stability or instability around a solution is of interest at early and also 
at late times. 

Then we will introduce also matter/radiation in the model and we shall generalise the method given in 
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H , Ref. [6], which permits to determine all critical points of a F{R) model. Ordinary matter is important in 

reconstructing the expansion history of the Universe and probing the phenomenological relevance of the 
models (see for example the recent papers [71 [HI H] , where the F{R) case has been discussed in detail). 
Our generalisation consists in the extension of that method in order to include all possible geometrical 
invariants. This means that F could be a generic scalar function of curvature, Ricci and Riemann tensors. 

To our knowledge, besides the paper |10j . the dynamical system analysis have been used and critical 
points derived and analysed mainly for models described by an action of the kind 

S = Sm + ^ f d^x^FiR), x = &t^Gn, (1.1) 



2X 

F{R) being an arbitrary function of the scalar curvature R, Gn the Newton's constant and 5*^ the matter 
action. The simplest choice F{R) = R — 2A corresponds to the ACDM model. Here we shall propose a 
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method which will permit to derive the critical points also for models described by a fmiction F[Uk), Uk 
being generic geometric (higher-order) invariants (see the Appendix). 

Their interest has been recently triggered by the appearance of Refs. [HHH] and after that investigated 
in many papers (see, for example [E]). Furthermore, very recently in [HI [HI [HI [ITl [18] , new "viable" 
models have been introduced and discussed. 

Also the stability of the solutions has been discussed in several places, an incomplete list being 
[ini[2ni[lIllllll3[Ml[ini[iail3IM[13 and earner references quoted in [Tg. To this aim, different 
techniques have been employed, including manifestly covariant and field theoretical approaches, where 
the gauge issue has been properly taken into account. All these investigations are in agreement with the 
following conditions which ensures the existence and the stability of the de Sitter solution with scalar 
curvature R — Rq: 



where F' and F" are the derivatives of F{R) with respect to R and everywhere in the paper the subscript 
indicates quantities evaluated at the critical point. The first condition in ()1.2p determines the scalar 
curvature of the de Sitter solution, while the second one gives the condition for the stability around such 
a solution. 

There are some theoretical (quantum effects and string-inspired) motivations in order to investigate 
gravitational models depending on higher-order invariants. The "string-inspired" scalar-Gauss-Bonnet 
gravity case F{R, G) has been suggested in Ref. [30] as a model for gravitational dark energy, while 
some time ago it has been proposed as a possible solution of the initial singularity problem f3T]. The 
investigation of different regimes of cosmic acceleration in such gravity models has been carried out in 
Refs. [3nilMI33l31[3Sl[Ml[37||3HllMllinillI]- in particular, in [35] a first attempt to the study of the 
stability of such kind of models has been carried out using an approach based on quantum field theory. 

The method we shall use in the present paper is based on a classical Lagrangian formalism [42, 43, 4^, 
inspired by the seminal paper [H] , where quantum gravitational effects were considered for the first time. 
With regard to this, it is well known that one-loop and two-loops quantum effects induce higher derivative 
gravitational terms in the effective gravitational Lagrangian. Instability due to quadratic terms have been 
investigated in [JS] . A particular case has been recently studied in [47i and general models depending on 
quadratic invariants have been investigated in [13 HH] ■ 

A stability analysis of nontrivial vacua in a general class of higher-derivative theories of gravitation 
has already been presented in [50]. Our approach is different from the one presented there since we are 
dealing with scalar quantities and moreover it is more general, since it is not restricted to the vacuum 
invariant submanifold. 

Finally, it should be stressed that the stability studied here is the one with respect to homogeneous 
perturbations. For the F{R) case, the stability criterion for homogeneous perturbations is equivalent to 
the inhomogeneous one [5Dj. 

The content of the paper is the following. In Section 2 we describe our general method in the vacuum 
case and derive general conditions for existence and stability of de Sitter solutions for arbitrary F{Uk) 
models. Then, in Section 3, we extend the method in order to include matter/radiation and give the 
conditions which (in principle) permit to derive all critical points. In Section 4 we analyse some specific 
models and finally we make some conclusions in Section 5. For the reader convenience, the paper ends 
with an Appendix in which we derive some interesting useful relations. 

2 The general case without matter 

In this Section we would like to study the stability condition of de Sitter solutions for models in which the 
Lagrangian density is an arbitrary function of all algebraic invariants built up with the Riemann tensor 
of the FRW space-time we are dealing with, that is 



where R is the scalar curvature, P — R^'^R^i, and Q = R""' ' Rap-jS are the two quadratic invariants and 
the dots means other independent algebraic invariants of higher order. 




2F(i?o) - RqF'[Rq) = , existence, 
^S^-1>0' stability. 



(1.2) 
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— F{R,P,Q,...), 



(2.1) 
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For the sake of convenience we write the metric in the form 

ds^ = -e2"(*)dt2 ^ e2a(t)^^2 ^ ^^^^ ^ g„(t) ^ ^^^^ ^ ga(t) (2.2) 

In this way Q!(i) ~ H[t) is the Hubble parameter and the generic invariant quantity U has the form 

U = e-2f "(*) M(n, d, a) = e-2p"(t) if, H) = H^p e-2p"(t) , (2.3) 

where 2p is the dimension (in mass) of the invariant under consideration and X = (H / — h/H) (see 
the Appendix). In particular one has 



R = 6e-2" 



2H^ + H -fiH 



(2.4) 



= 6i?^e-^"(2 + X), 

P = 12e-4" Ih'^H^ - 3ni7^ - 2nHH + 37?^ + 'iH'^H + ij^") = 12ij4e-4"(3 + 3X + X^) 
Q = 12e-4" hi^H'^ - 2hH^ - 2niJij + 2H'^ + 2H'^H + hA = UH*e-'^''{2 + 2X + X^) , 
G^Bj^ -AP + Q = 24e~''"(i?'' + H'^H) = 2m^e^'^''{l + X) . 

Using this notation, the action reads 

S^- I (fix [ dtL{n,h,a,d,a) ^ — ( d^x [ dt e"+^" F{n,n, a, a) (2.5) 



2X 

and the Lagrange equations corresponding to the two Lagrangian variables n{t) and a(t) are given by 

i^„ = |^-|^ = 0, (2.6) 
on at on 

+ (2.7) 

" da dt da dt'^ da 

Since the function n{t) is a "gauge function" , which corresponds to the choice of the evolution parameter, 
we expect equation (|2.7p to be a direct consequence of equation (j2.6p and in fact a straightforward 
calculation leads to the identity (see Eq. (jA.Sp in the Appendix) 

^ = nE,, + aE^ . (2.8) 
dt 

As a consequence of such an identity we may limit our analysis to the use of equation (|2.6p , much simpler 
than equation (j2.7p . Furthermore, we may use the gauge freedom and fix the cosmological time by means 
of the condition N{t) = 1, that is n{t) — 0. From now on it is understood that all quantities will be 
evaluated in such a gauge and so the parameter t corresponds to the standard cosmological time. In this 
way Eqs. ([^ and (P?7)) read 

F + Fn- 3HFn HFjj - HFh + F^ HF^ + 3H^Fj^ = , (2.9) 

F^^Fh+ 6HFfj - 3HFh + 3F + iHFfj + QH'^Ffj ^ , (2.10) 



where Fx — dxF means derivative of F with respect to X. All derivatives with respect to n in Eq. (j2.9p 
have been eliminated by using Eq. (jA.6p in the Appendix. 

In the next Section we shall study all critical points of the system in the presence of matter, while 
here we shall limit ourselves to determine the conditions for the existence and the stability of de Sitter 
solutions in the absence of matter. To this aim we take into account Eq. ()2.9p and re-write it as a 
differential equation for the Hubble parameter H{t), that is 

HF„^ H + H{F„^ -Fh)H + S(F, H) = , (2.11) 

where we have set 

E{H,H) = F - HFh + 'iH^Ffj . (2.12) 
We see that the condition for the existence of a de Sitter solution Pq = {H — 0, H — Hq) is 

E(i7o,0)=0 =^ [F-HFh + 3H^F^]p^^^0, (2.13) 
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where here all quantities have to be evaluated at Pq. The latter condition, which gives rise to the de 
Sitter critical point, can be simplified by means of Eq. (|A.9p in the Appendix. In fact we have 



[F - HFh + SH'Ff,] = [F- H^Ffj] = . 



(2.14) 



In order to study the stability of such a solution we have to distinguish between two possible cases. The 
simplest one occurs when F is linear in H , that is Fjjjj = 0. In such a case Eq. ()2.1ip assumes the form 



{HF^^ - HF^ - + 0) = 0) 



H = Ho. 



(2.15) 



Then we see that in such a special case the solution, if it exists, is trivially stable, since the field equation 
reduces to an algebraic equation which fixes the value of H in terms of the parameters of the system. 
The Einstein's theory with cosmological constant = i? — 2A is an example of this kind and in fact 
Eq. ((TT5)) gives rise to H"^ = = A/3. 

In the second case F^^ ^ 0, in order to study the stability one may transform Eq. (|2.1ip in an 
autonomous system by introducing the functions K(t) = H{t). In this way, one gets 



H = K, 



K 



m,K) + iHF^^-Fjj)K 



HF, 



(2.16) 



HH 



The critical points are determined by the condition (H = 0, K = 0) and the linearized system which 
determines the stability reads 



\5K 



M = 



1 

^0 Bo 



where 



Ao = 



HF, 



HH 



Fhh — iHF 



HH 



6F, 



H 



HH 



(2.17) 



(2.18) 



Po 



Bo 



Ph - H^HH - ^H 



HF, 



HH 



-Ho. 



(2.19) 



The stability conditions are obtained by requiring TrM < and detAf > 0. The first condition is 
trivially satisfied since Ho is positive, while the second one gives 



3HF„rr - Fhh + 6F, 



H 



F, 



HH 



>0, 



Summarizing, we have the general results 

F — H'^Ffj] = , critical points, 

> , stability condition when F^ ^ 



3HF„^-Fhh+6F^ 



^0, 



(2.20) 



(2.21) 



In the particular case in which F^^ = the eventual (de Sitter) critical point is always stable. Of course, 
when F — F{R), Eqs. (|2.2ip become equivalent to Eqs. ([TT 



3 The general case with matter 

The inclusion of matter in the model is obtained by adding the matter action Sm to the gravitational 
action (j2.5p . Of course, in order to preserve symmetry we have to consider homogeneous and isotropic 
matter, that is a perfect fluid. First of all we observe that 

^^'" 6n = - [ d^x^gToog'^Sn, (3.1) 



5n 
6S„ 



Sa = - (Tx^/^Tabg 5a, (3.2) 



5a 

where T^^ = (Too, Tab) is the energy-momentum tensor, which for a perfect fluid satisfies the conditions 
W° = -P, Tabg^' = 3p, a, 6 =1,2, 3, 
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p and p — p{p) being respectively the density and the pressure of matter. 

Now Eqs. (|2.6p and (|2.7p triviaUy changes for the presence of matter and in fact we get 



^ _dL d dL 
" dn dt dri 



Ea — 



dL d dL 



^V-gToog 
rf2 dL 



00 



da dt da dt^ da 



= \r-9Tab9 



ab 



'6p 1 



'9- 



The identity (|2.8p is stiU vahd and is equivalent to the energy-momentum conservation law 

V^T^'' = =^ p=-3H{p + p). 
Using the notation of the previous Section and putting again n(t) = now we get 

HFj^ - HFh +F- HFfj + ■iH'^Ffj = 2p , 



(3.3) 
(3.4) 

(3.5) 
(3.6) 



F^^Fh+ 6HFfj - 3HF„ + 3F + iHF^ + %H'^F^ 



. 



(3.7) 



The latter equations are the generalisation to arbitrary action of the well known Friedmann equations. 
It is interesting to observe that in the pure Einstein gravity, that is for F = R, they read 



H^Ffj ^ F: 



X 



2p 



1. 



{'3H^ + 2H)Ffj = (3 + 2X)Fx 



-6p 
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where we have introduced the dimcnsionless variables 
2p 2p_ 2p _2p 

Fx ' " " H^F„ Fx 



~ H^F 



H 



(3.8) 
(3.9) 

(3.10) 



which in this special case are given by the usual values flp = p/3H^ and fip — p/3H^. From Eqs. 
and (|3.9p it follows 



p lip 6 



(3.11) 



In the general case, Eqs. p.6p and (|3.7p have more terms with respect to p.Sp and p.9p and it is quite 
natural to interpret them as corrections due to the presence of higher-order terms in the action. Then 
we define 



^curv 



1 



^curv 



-1 



3 



(3.12) 



Weff 



eff 



3 



(3.13) 



where ri™*"" and Jl™'"" are complicated expressions, which only depend on the function F. They can be 
derived from (j3.6p and p.7p . but they explicit form is not necessary for our aims. The effective quantity 
Wef f is equal to the ratio between the effective density and the effective pressure and it could be negative 
even if one considers only ordinary matter. It is known that the current-measured value of Weff is near 
-1. 

In order to get all critical points of the system now we follow the method described in [5j . First of all 
we introduce the dimcnsionless variables 



H 



2p_ 

Fx 



H 



2p_ 

Fx 



(3.14) 



X 



H 



Y 



F ~ HFh _ F 
~lPF~ " F^ 



X . 



Z 



HF 



H 



Fx 



(3.15) 
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where the prime means derivative with respect to a and the quantity 

e = = = (3.16) 

has to be considered as a function of the variables X and Y. In general it is a function of X and H, but 
this latter quantity can be expressed in terms of X and y as a direct consequence of the definition of Y 
itself. Then we derive an autonomous system by taking the derivatives of such variables. From Eq. (|3.6p 
we have the constraint 

np = Y + Z + 3 =^ ncurv = -{Y + Z + 2), (3.17) 

which reduces to the standard one when F is linear in R (general relativity with cosmological constant). 

Deriving the variables above by taking into account of (|3.5p and p.7p we get the system of first order 
differential equations 

X' = -2X^ - jX + f3(Z + £,) , 
Y' = -{2X + Z + C)Y - XZ. 

Z' = -Sflp ~ [X{Z + 6) + Z(Z + e + 6) + 3(r + ^ + 3)] , (3.18) 

n'p = -3(f7p + -{2X + Z + c)np , 

O; = -3^ {Qp + Qp) -{2X + Z + , 

where X' = ^ = ^'^ -"^^ ^^"^ P ~ P^P^ been assumed. The first equation in (|3.18p has been 

obtained by deriving the function F^ with respect to t and putting 

F F F J-f F 

P^P{X,Y)^-^^-^, ^^^^x,Y)^-f^^—^^P^x+i. (3.19) 

It is understood that F^^ ^ has been assumed. 

Note that in principle one could consider a mixture of different kind of matter/radiation with den- 
sities pk and corresponding pressures pk- In such a case in Eqs. p.lSp there will appear the sums of 
the corresponding quantities fip^. and Op^. and the corresponding differential equations for any type of 
matter/radiation. For simplicity here we consider only one type of matter/radiation described by the 
equation of state p = p{p) . 

In general, among the five differential equations (j3.18p only three of them are linear independent for 
the presence of the constraint (j3.17p and the equation of state. For simplicity now we assume the pressure 
to be proportional to the density, that is p = wp, with constant w, where for ordinary matter < w < 1/3 
{w = corresponds to dust, while w = 1/3 to pure radiation), but in principle one could also consider 
"exotic" matter with w < and cosmological constant which corresponds to w = —1. With this choice 
and from Eq. (|3.17p we get 

rip = wnp , np = z + Y + 3. (3.20) 

As a consequence the autonomous system which gives rise to the critical points can be chosen as 
= X' = -2X^ -"fX + P{Z + 

= Y' = -{2X + Z + ^)Y - XZ (3.21) 
= Z' = -3(1 + w){Z + Y + 3)-{Z + C){Z + 3)- X{Z + 6) 

and rip at the critical points will be determined by means of Eq. ()3.17p . The critical points are the 
solutions of the algebraic system p.2ip . The number and the position of such points depends on the 
Lagrangian throughout the functions /?, 7 and ^. In principle, given F one can derive all critical points 
as the constant solutions of Eqs. (|3.2ip . but in practice for a generic F the algebraic system could be 
very complicated and the solutions quite involved. We shall consider in detail some particular cases in 
the next Section. 

A brief comment about the equivalence of (|3.2ip and Eqs. p.6p - p.7p is in order. It is well known that 
the Gauss-Bonnet scalar G does not contribute to the field equations, because it is a topological invariant 
and this means that the field equations (|3.6p - p.7p are invariant with respect to the transformation 
F ^ F + cG, c being an arbitrary constant. On the contrary, the new variables X, y, Z (by definitions) 
are not invariant with respect to such a transformation and so the system p.2ip could contain the 
parameter 7 explicitly, but nevertheless the solutions for H{t) do not depend on such a parameter. The 
same thing is true for the critical points. 
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It has also to be noted that in general, given a Lagrangian F, there are solutions of the system (|3.2ip . 
which give rise to trivial or unphysical models. Of course such solutions have to be dropped (see specific 
examples) . 

Before to analyse all possible solutions of ()3.2ip . we observe that in the absence of matter the de 
Sitter solution X = oi course is a critical point for the system if the first equation in ()2.21|) is satisfied. 
In fact we trivially see that X = 0, ^Ip = Y + Z + 3 = satisfy (|3.2ip if F — ^ + 3 = 0, which is equivalent 
to first equation in (|2.2ip when X ~ 0. (Strictly speaking, the condition F = 1 is equivalent to (|2.2ip . 
when Hq ^ 0. In principle, Eq. (|2.2ip can also have the Minkowskian solution i?o = 0). On the de Sitter 
solution the value ofweff is exactly equal to —1. By taking the variations of p.2ip with respect to X and 
Y on the critical point and using the identity (|A.9p in the Appendix we also get the stability condition 

in arm . 

In order to show that the present approach is a generalisation of the one proposed in [Sj, we observe 
that when the function F ~ F{R) depends on the unique invariant U — R, we exactly obtain all known 
results. For such a special case, using ()2.4p (with n{t) = 0) and (|A.5P - (jA.6p in the appendix we get 

^ HF,^2pU 

^ Fx Ux ^ ' 

^ " " FuUxx + FuuUl " (SmFRR ' ^^'^^^ 

7 = ^^?^=/3Cx+e = 4, (3.24) 

Now with the trivial replacement 

X = X3-2, y==a;2-X3 + 2, Z = -(xi+4), = — , (3.25) 

m 

Eqs. p.lSp become equivalent to the analog equations in Ref. H], where this class of models have been 
studied in detail. 

As already said above, the critical points of (|3.2ip depends on ^, /3, 7, which in general are complicated 
functions of X and Y , then it is not possible to determine general solutions without to choose the model, 
nevertheless it is convenient to distinguish two distinct classes of solutions characterised by the values of 
w ^ —1 and w = — 1. For the sake of completeness we consider w < 1/3 and so we write the solutions 
also for "exotic" matter, that is quintessence (—1 < w < 0) and phantom {w < —1). Of course, such 
solutions have to be dropped if one is only interested in ordinary matter/radiation. We have 

• w ^ ~1 — The critical points are the solutions of the system of three equations 

2X^ +jX-(3{Z + = 

{2X + Z + S,)Y + XZ = (3.26) 
3(1 + w){Z + r + 3) + (Z + e)(^ + 3) + X{Z + 6) = 

where ^, /3, 7 are functions of X, Y determined by Eqs. (|3.16p and p.l9p . The stability matrix has 
three eigenvalues and the point is stable if the real parts of all of them are negative. 

The latter system has always the de Sitter solution Pq = {X = Q,Y = 1, Z = —4), where flp = 
and Weff = — 1. Note however that such a solution could exist also in the presence of matter, since 
the existence of Pq critical point only implies that the critical value for flp vanishes. 

51p^0,w = — 1 — The critical points are given by 

2X^ +jX-l3{Z + = 

{2X + Z + ^)Y + XZ = w = -l. (3.27) 

{Z + C){Z + 3)+X{Z + 6) = 

For this class of solutions, the non-singular stability matrix has three eigenvalues and the point is 
stable if the real parts of all of them are negative. 

We see that there is at least one singular case (critical line) when X = and Z = — ^ = —4. In fact 
in such a case Y or f2p are undetermined since 

= F + 3-C(0,y) = F- 1 =^ Y = l + np, arbitrary. (3.28) 
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Such a solution can be seen as a generalisation of the de Sitter solution for a model with cosmological 
constant. The de Sitter critical point for the model F = F — 2K reads {X = 0, F = l^Z = —4). 
Such a solution follows from Eq. (13.28|) if we choose po = A. In fact, on the critical point {X = 
0, r = 1 + r^p, Z = -4) (Eq. (I32H1)) and from definitions ([3l^ we get 

"In F _ P 
n„ = -^f— = — 1 =^ Y= ^ = 1 , (3.29) 

which corresponds to de Sitter critical point for F. Of course, Eq. (j3.28p is more general than the 
case with pure cosmological constant since p is not necessary a constant. 

Of course for this special class of solutions w^^jf = — 1. Note also that the stability matrix has always 
a vanishing eigenvalue and the stability of the system is determined by the other two eigenvalues. 

For some models, but just for technical reasons, it could be convenient to treat the cosmological 
constant as matter, using the previous identification we have done. 



4 Explicit examples 

In order to see how the method works, now we give explicit solutions for some models and, when possible 
we also study the stability of the critical points. We restrict our analysis to the values < w < 1/3 and 
to the special value w = — 1, which corresponds to the pure cosmological constant, but in principle any 
negative value of w could be considered, even if this will be in contrast with the aim of modified gravity. 
In fact, modified gravity can generate an effective negative value of w without the use of phantom or 
quintessence. 

It as to be stressed that in general, due to technical difficulties, one has to study the models by a 
numerical analysis. Only for some special cases one is able to find analytical results. Here we report the 
results for some models of the latter class in which the analytical analysis can be completely carried out. 
We also study more complicated models and for those we limit our analysis to the de Sitter solutions. 

In the following we shall use the compact notation 

P={X,Y,Z,np,w,ff), Po = (0, 1,-4, 0,-1), PA-(0,l + r!A,-4,r!A,-l). (4.1) 

The latter is an additional critical point that we have for the choice w — —1 and can be seen as the de 
Sitter solution in the presence of cosmological constant. 

F = R — fi'^/R — This is the well known model introduced in [TTl [T^] and discussed in [S]. For this 
model the system p.2ip with arbitrary w has six different solutions, but only two of them effectively 
correspond to physical critical points, if < w < 1/3. In principle there are other critical points 
for negative values of w (phantom or quintessence) and moreover there is also a particular solution 
for w ~ —I which corresponds to the model with a cosmological constant A. 

Solving the autonomous system one finds 

• P = Pq: unstable de Sitter critical point. The critical value for the scalar curvature reads 
Ro = V3Ai'. 

• P = (— 1/2, — 1, — 2, 0, —2/3): stable critical point. At the critical value, iJo ~ 0. 

• P = (3(1 + w)/2, -(5 + 3w), -2(5 + 3ti;), -3(4 + 3ti;), -(2 + w)): unstable critical point where 
Ho =0. 

• P = P\: unstable critical point. At the critical value one has Hq — (A/6)(l + -y/l + 3/i^/4A2 . 

F = R + aR? + bP + cQ — (Starobinsky-like model). Here we have to assume 3a + b + c ^ Q otherwise 
the quadratic term becomes proportional to the Gauss-Bonnet invariant. For < w < 1/3, this 
model has only one critical point. In order to have a de Sitter solution, we have to introduce a 
cosmological constant A. We have in fact 

• P = Pf)-. Minkowskian solution with _Ro = 0, which is stable if 3a + 6 + c > 0. 

• P = Pa'- de Sitter critical point with _Ro = 6A, which is stable if 3a + + c > 0, in agreement 
with [ig. 
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F = R — (PQs — This is the simplest toy model with the cubic invariant = R'^^''^ RaPfiuRtl^s ■ For 
this model we have the following critical points: 

• P = Pq: unstable de Sitter solution with Rq = 6/d. 

• P ~ Pq: stable Minkowskian solution with Rq = 0. 

• P = (0.05, 0.60, -3.60, 0, -1.03): stable solution with Hq = 0. 

• P = P\: this point exists, but not for any value of d and A. Also the value of Hq and the 
stability depend on the parameters. 

F = R + aR' + bP + cQ ~ d'Q^ — This is a generalisation of the previous two models. It may be 
motivated by the two- loop corrections in quantum gravity |5H I52j . The de Sitter critical points 
have been studied in Ref. [3j. The algebraic equations p.21|) are too complicated to be solved 
analytically, but it is easy to verify that there are at least the following solutions: 

• P ~ Pq: de Sitter solution with Rq — 6/d. This is stable if 3a + 6 + c + 3(i > 0. 

• P ~ Pq: Minkowskian solution with Rq = 0, which is stable if 3a + 6 + c > 0. 

• P — Pa: also in this case this point exists and is stable depending on the parameters (see 



5 Conclusion 

In this paper we have presented a general technique which permits to arrive at a first order autonomous 
system of differential equations classically equivalent to the equations of motion for models of modified 
gravity based on an arbitrary function F(R, P,Q,Q3...), namely built up with all possible geometric 
invariant quantities of the FRW space-time. Dynamical system techniques have been applied to the 
investigation of critical points. We have shown that, in the special case of F{R) theories, the method 
gives rise to the well known results O [71 [8] , but in principle it can be applied to the study of much more 
general cases. 

As applications, we have considered some simple models, for which a complete analytical analysis 
have been carried out. However, in general, due to technical difficulties, a numerical analysis is required. 
Among the models investigated, we would like to remind that we were able to deal with one which involves 
a cubic invariant in the curvature tensor and, to our knowledge, this has never been considered before, 
and this shows the power of our approach. 



A Appendix 

Here we show that in the cases we are considering, Eq. ()2.8p is always satisfied. To this aim, as in Eq. 
we denote by U{x) a generic invariant quantity of dimension 2p (in mass) computed in the system of 
coordinates {x} = {{t,x)} and by U{y) the same quantity computed in the system {y} = {(t, y)}. Such 
coordinates are chosen in such a way that dr = e"(*)dt and as a consequence the metrics read 



ds^ 



(A.l) 



Recalling that U is built up with the Riemann tensor, for dimensional reasons in the system {y} one has 



Uiy) = E 



k=0 



da{T) 




'd'air)' 






dr' 



and since under diffeomorphism [/ is a scalar quantity we also have 



k=Q 



'da{t)' 


2{p-k) - 


dt 





d'a{t) da{t) dn{t) 
dt^ Jt df 



n k 



(A.2) 



(A.3) 



This means that in a generic system of coordinates an arbitrary invariant can be written in the form 

a n H 

AfcA-, X = 

k=0 



a 



n 



(A.4) 
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From the latter equation we directly get 



TJ. - _Ux. 

u- - ^ 



Uh 



-2pU, 

2pU _ 2HUx 
H m 



n Ux 



and for a generic function F{Ui, U2, U3, ...Ua, •■■) 

From equations above it directly follows 
F„ + HFh + 2HFjj +nFn^{). 



p. - _£x 

p. - Fx 



2XFx 
H 



fiFx 
IP' 



(A.5) 



(A.6) 



(A.7) 



After these considerations it is quite easy to recover Eqs. p.8|l . In fact, using Eqs. p.6p . p.7p . (|A.6P - (|A.7|) 
we get 



En - flEn - aEa = 



[fi + 3a){Fn + HFh + ^HF^ + riFn) 



HFh + 'iHFfj + iiFn 



Before to end the appendix we also derive the useful relation 

(X=0,n=0) 



HE 



H 



HOh log ^ 



1. 



(X=0,n=0) 



The latter identity can be derived as follows. From (|A.6|) we have 



Eh 
HE 



H 



HFh 



Ex 



{X=0,n=Q) 



(A.; 



(A.9) 



(A.IO) 



where Ak{Ua) are the coefficients of the invariant Ua as in (jA.4[) . Then eq. (|A.9p is true if for any invariant 
Ua one has 



2A^{Ua)^PaAo{Ua). 



(A.ll) 



Such a relation is a direct consequence of the form of the Riemann tensor. In fact, in this case the non 
vanishing components read 



Ra 



hah 



a, 6 = 1,2,3 



and since any invariant is built up with Riemann tensor it is easy to sec that the relation (jA.lip is always 
satisfied. 
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